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We propose a new combined approach of the exact diagonalization, the renormalization group 
and the Bethe ansatz for precise estimates of the charge gap A in the one- dimensional extended 
Hubbard model with the onsite and the nearest-neighbor interactions U and V at quarter filling. 
This approach enables us to obtain the absolute value of A including the prefactor without ambiguity 
even in the critical regime of the metal-insulator transition (MIT) where A is exponentially small, 
beyond usual renormalization group methods and/or finite size scaling approaches. The detailed 
results of A down to of order of 10~ 10 near the MIT are shown as contour lines on the U-V plane. 

PACS: 71.10.Fd, 71.27.+a, 71.30.H-h 



There has been much theoretical interest in the one- 
dimensional (ID) strongly correlated electron systems 
such as the t-J model and the Hubbard model as a 
good testing ground for the concept of the Tomonaga- 
Luttinger liquid 1,2 . Various methods, such as the weak 
coupling theory (g-ology), the bosonization theory, the 
Bethe ansatz (BA) method, the conformal field theory 
and the numerical approaches have been used to clarify 
the nature of these models 3 . Among them, combined 
approaches of the exact diagonalization (ED) and the 
renormalization group (RG) methods have been exten- 
sively developed to investigate the critical behavior near 
the quantum critical point 4-9 . These approaches enable 
us to obtain accurate results of the phase boundary for 
the spin gap phase and those for the charge gap phase 
beyond the purely numerical approaches combined with 
the usual finite size scaling. 

Recently, we have intensively examined the critical be- 
havior near the metal-insulator transition (MIT) in the 
one-dimensional extended Hubbard model with the on- 
site and the nearest-neighbor interactions U and V at 
quarter filling using a combined approach of the ED 
and RG methods 7-9 . In this approach, the Luttinger- 
liquid parameter K p is calculated by using the ED for 
finite size systems and is substituted into the RG equa- 
tion as an initial condition to obtain K p in the infinite 
size system. The obtained result agrees very well with 
the available exact result for U = oo even in the crit- 
ical regime of the MIT where the characteristic energy 
becomes exponentially small and the usual finite size 
scaling is not applicable. When the system approaches 
the MIT critical point V -> V c for a fixed U, K p be- 
haves as (K p — i) 2 = Cif(l — V/V c ), where the criti- 
cal value V c and the coefficient ck are functions of U 9 . 
This approach also yields the critical behavior of the 
charge gap A in the insulating state near the MIT, where 
| In A| -2 = ca(V/V c — 1) with the coefficient ca which is 
a function of U. In these studies 7-9 , however, the ab- 



solute value of A including prefactor was not explicitly 
obtained. 

In general, it is considered to be difficult for the RG 
method or its derivative method to yield absolute value 
of physical quantities including prefactor. To overcome 
this difficulty, we join the BA method to our previous 
combined approach of the ED and the RG methods in 
the present study. More explicitly, the BA result in the 
infinite size system with U = oo is connected to the ED 
result in the finite size system with finite U through the 
analysis of the RG solution. The new combined approach 
enables us to estimate the absolute value of A including 
the prefactor without ambiguity in contrast to the previ- 
ous combined approaches. 

The extended Hubbard model is given by the following 
Hamiltonian 



H = -tJ2i<&rCi+i* + h.c) 



(1) 



where cL stands for the creation operator for an elec- 



tron with spin a at site i and rij 



t repre- 



sents the transfer energy between the nearest-neighbor 
sites and is set to be unity (t=l) in the present study. 
It is well known that this Hamiltonian eq. (1) can be 
mapped on an XXZ quantum spin Hamiltonian in the 
limit U — > oo. The term of the nearest- neighbor inter- 
action V corresponds to the Z-component of the anti- 
ferromagnetic exchange coupling and the transfer energy 
t corresponds to the X-component of that. When the 
Z-component is larger than the X-component, the sys- 
tem has a "/sing" -like symmetry and an excitation gap 
exists. For the Hubbard model, this corresponds to the 
case with V > 2t where the exact result of the charge 
gap is given by 10 
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A = 4(sinhA) J2 
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2 coshnA 



with 



A = - In(V/2 - v/(F/2) 2 - 1)). 



(2) 



(3) 



On the other hand, in the case with "XF"-like symmetry 
(V < 2i), the system is metallic and the Luttinger-liquid 
parameter K p is exactly obtained by cos(jg-) = — V/2 11 . 

In order to introduce our approach, we briefly discuss 
a general argument for ID-electron systems based on the 
bosonization theory 1-3 . According to this theory, the 
effective Hamiltonian for the ID electron systems can 
be generally separated into the charge and spin parts. 
Therefore, we turn our attention to only the charge part 
and do not consider the spin part in this work. In the 
low energy limit, the effective Hamiltonian of the charge 
part is given by 



£ dxcoB[2y/84> p (x)], (4) 



+ 



where v p and K p are the charge velocity and the coupling 
parameter, respectively. The operator (j) p and the dual 
operator 9 p represent the phase fields of the charge part. 
<?3_l denotes the amplitude of the umklapp scattering and 
a is a short-distance cutoff. On the basis of the Hamilto- 
nian eq. (4) , the electronic state is described by only the 
two parameters K p and g^± except for the energy scale 
determined by v p . 

At quarter filling, the 8kp umklapp scattering plays 
the crucial effect for the charge gap. The effect of the 
umklapp term is renormalized under the change of the 
cutoff a — > e e a, where £ is the scaling quantity. This 
process is also considered as the change of the system 
size L — > e e L. Therefore, the size dependence of K p is 
described by the RG equations. 7-9 In this work, we adopt 
the Kehrein's formulation as the RG equations 13 ' 14 



dK p {£) 

d£ 

dlogG(l) 
dt 



G{£fK p {£f 

T{&K p {t)-iy 

= [2 -8K p (£)}, 



(5) 
(6) 



where T(x) is T-function and G(l) stands the umklapp 
effect with G(0) = g3±/(2wv p ). Here, the value of the 
short-distance cutoff a is selected to a lattice constant 
of the system and set to be unity. This formulation is 
an extension of the perturbative RG theory and allows 
us to estimate the charge gap beyond the weak coupling 
regime. 

To solve the RG equations concretely, we need an ini- 
tial condition for the two values: K p (0) and G(0). Be- 
cause it is easy for the ED calculation to obtain K p {t) 



as compared to G(£), we eliminate G(£) in the RG equa- 
tions. For this purpose, we integrate eq. (6) to yield 



G{l) = G{i 1 )J>- &K ^' )]M \ 



(7) 



where l\ is a constant. Substituting eq. (7) into eq. (5), 
we obtain the differential equation for K p (£) as 



dK p (l) _ GHt^J^^K, 



{If 



de 



T(8K P (£) 1) 



(8) 



Setting K p (ii) as the initial condition, we solve eq. (8) 
numerically except the constant G{£\). The value of 
G(£i) is determined by comparing the solution K p {£) at 
£ = I2 with the initial value K p (£ 2 )- Then, the solutions 
for eqs. (5) and (6) are completely obtained. 

Using the relation £ ~ hiL, we calculate two initial 
values K p (£i) and K p (£ 2 ) with L\- and L2-site systems 
by the ED method. In the finite size systems, K p (L) is 
calculated by the charge susceptibility \c and the Drude 
weight D by 



K p = ±(n Xc D) 1 / 2 



(9) 



with D 



n d 2 E„(4>) 



T — dip — ' wri ere Eq((/)) is the total energy 
of the ground state as a function of a magnetic flux Lcj) 
and L is the system size 3 . Here, the magnetic flux is im- 
posed by introducing the following gauge transformation: 
for an arbitrary site m. The uniform 



4 



charge susceptibility \c is obtained from 

4/L 

Xc ~ E (N + 2,L) + E (N - 2, L) - 2£ (A, T) 



(10) 



where Eq(N,L) is the ground state energy of a system 
with L sites and N electrons. Here, the filling n is defined 
by n = N/L. We numerically diagonalize the Hamilto- 
nian eq. (1) up to 16 sites system using the standard 
Lanczos algorithm. In the case with U = 00, we also 
calculate Eo(N,L) by using the Bethe ansatz method 12 
for finite size systems up to 800 sites system. Using the 
definitions of eqs. (9) and (10), we calculate K p and \c 
from the ground state energy of the finite size system. 

In Fig. 1, we show the size dependence of K p (L) ob- 
tained from the solutions of the RG equations together 
with the exact Bethe ansatz results for various V at 
U = 00. We choose L\ = 12 and L 2 = 16 for the nu- 
merical initial condition in the RG equations, and we set 
Li = L 2 - 4 hereafter. The limit K p (L — > 00) of the BA 
result becomes a finite value for V < 2 and converges to 
zero for V > 2. We see that the RG solution is very close 
to the exact result and the size dependence of K p (L) is 
well described by the RG equations except very large size 
systems for V > 2. 

In Fig. 2, we show the RG flows of the systems with 
various V for U = 10 together with those for U = 00. 
We observe that the RG flows for U = 10 and those for 
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FIG. 1. The size dependence of K P (L) obtained from the 
RG equations (solid lines) and the exact Bethe ansatz results 
(filled circles) for various V at U = oo. The shadowed squares 
are the numerical initial conditions with L\ = 12 and L2 = 16 
for the RG equations. 

U = 00 coincide to each other. Based on the Luttinger 
liquid theory, it is considered that the systems assigned 
by the same RG flows have the same electronic state ex- 
cept energy scale v p . If we find out a RG flow for U = 00 
(we call it a reference system) corresponding to that for 
finite U (it is a original system), we can connect both sys- 
tems through the RG flow and derive properties of the 
original system from the known result of the reference 
system. To identify the RG flow, we use a renormal- 
ized coupling constant G{1) constructed by the product 
of G(£) and the effective energy scale e~ e ( 2 ~ 8Kp \ In the 
limit £ — > 00, G(tj diverges in proportional to e £ ( 2 ~ 8Ap ) 
(see eq. (7)), while G(oo) remains a finite value and it 
becomes a unique index characterizing the RG flow. 

Here, we determine the nearest neighbor repulsion V 
in the reference system with U = 00 so as to fit the RG 
flow of the reference system to that of the original system 
with original's V for a finite value of U. The reference's 
parameters V corresponding to the several original's pa- 
rameters V for U = 10 are shown in Table I. The ref- 
erence's V is smaller than the corresponding original's 
V. This suggests that the on-site repulsion U causes 
the renormalization of the nearest neighbor repulsion V. 
Fig. 2 also shows that the point indicating the initial con- 
dition of the reference system is located at downstream 
than that of the original system on the RG flow. This 
means that the effective size of the reference system is 
larger than the original system size. 

Substituting the reference's V into eq. (3), we obtain 
the charge gap A r of the reference system from eq. (2). 
Taking into account the difference of the energy scale, 
i.e., the charge velocity between the original and the ref- 
erence systems, we estimate the charge gap of the original 




at U = 10 with the numerical initial condition (filled cir- 
cles) . The open squares are the numerical initial condition at 
U = 00 fitting to the RG flow for U = 10. 



system as A = v°/v p A r , where v° and v p are the charge 
velocities of the original and the reference systems, re- 
spectively. Both of the charge velocities are calculated 
from the ED results with same size systems through the 
relation of v p = D/2K p . 

In Fig. 3, we plot the charge gap A as a function of 
V - V c at U = 10 for L 2 = 8, 12 and 16 together with 
the result of an extrapolation with L 2 — > °°- Here, V c 
is the critical value of the MIT and determined by the 
condition K p (l — > 00) = 1/4. The values of V c are given 
by 2.684, 2.583 and 2.567 for finite size systems with 
L2 =8, 12 and 16, respectively, which yield an L2 — ► 00 
extrapolated value V c — 2.55. We note that detailed 
analyses of V c and the MIT have been already discussed 
in the previous works 6-8,15-17 . The -^-dependence of A 
is assumed to be proportional to I/L2, resulting in an 
extrapolated value of A with Li — > 00 as shown in Fig. 
3. The obtained result of A is in good agreement with 
the recent DMRG result of A at V - V c ~ 1.45 18 . The 
inset in Fig. 3 shows the semi-log plot of A near the MIT. 
The system size dependence of A is very small even in the 
critical regime of the MIT near V c . These results show 
that the new combined approach is especially efficient to 
analyze the very small charge gap near the MIT. 

In Fig. 4, we show the detailed results of contour lines 



TABLE I. The original's and the reference's parameters V 
for U — 10 (see in the text). 



Original's V 


1.0 


1.5 


2.0 2.5 


3.0 


3.5 


4.0 


Reference's V 


0.38 


0.85 


1.33 1.84 


2.38 


2.92 


3.43 
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FIG. 3. The charge gap A as a function of V- V c at U = 10 
for L 2 = 8, 12 and 16 together with the extrapolated result 
(dashed line). The solid diamond is the DMRG result 18 . Inset 
shows the semi-log plot of A near the MIT. 
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FIG. 4. The contour lines of the charge gap A near the 
MIT on the U-V plane. The solid line represents the phase 
boundary of the MIT. 



for the extrapolated value of the charge gap A near the 
MIT on the U-V plane. Rough estimation of A has been 
already reported in our previous paper 9 . However, the 
previous result was limited in the case with large gap 
(A > 0.25) region, since the usual finite size scaling is 
used for the ED result. We stress that our new combined 
approach has an ability to estimate the very small charge 
gap of order of 10 -10 far from the limitation of the usual 
numerical estimation. 

In summary, we examined the new combined approach 
of the ED, the RG and BA methods to clarify the charge 
gap A of the ID extended Hubbard model with the on- 
site and the nearest-neighbor interactions U and V at 
quarter filling. Analyzing the solution of the RG equa- 
tions, we connect the original system with a finite U to 
the reference system with U = oo in which the charge 
gap has obtained as a function of V by using the exact 
BA method. Adjusting the parameter V of the reference 
system so as to fit the RG flow of the reference system 
to that of the original system, we estimate the absolute 
value of A of the original system including the prefactor. 
This approach is able to supply us with unambiguous 
and accurate result of A beyond the usual RG method 
and/or the ED method, even if energy scale becomes ex- 
ponentially small. Detailed Analysis of A is shown as 
the contour lines on the U-V plane in the critical regime 
near the MIT with very small gap. 
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